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ABSTRACT 


Observed  ocean  acoustic  wavefronts  show  surprising  stability  in  long-range 
acoustic  transmission  experiments.  This  suggests  that  ocean  scattering  processes  tend  to 
redistribute  acoustic  energy  dominantly  along  the  wavefront  rather  than  across  it.  The 
purpose  of  this  thesis  is  to  elucidate  the  physical  mechanism  for  this  type  of  scattering  by 
presenting  a  ray-based  physical  model.  Analytic  fonnulae  are  presented  that  predict 
wavefront  distortions  caused  by  ocean  internal  waves  and  other  processes.  Applications 
of  this  study  include  wavefront  healing  near  underwater  obstacles,  out-of-plane  scattering 
and  the  vertical  redistribution  of  energy  of  off-axis  sources. 
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I.  INTRODUCTION 


Observed  ocean  acoustic  wavefronts  show  surprising  stability  in  long-range 
acoustic  transmission  experiments.  This  suggests  that  ocean  scattering  processes  tend  to 
redistribute  acoustic  energy  dominantly  along  the  wavefront  rather  than  across  it.  The 
physical  mechanism  for  this  type  of  scattering  has  not  yet  been  elucidated. 

There  are  several  practical  implications  for  this  type  of  scattering.  First, 
scattering  along  the  wavefront  leads  to  the  important  result  that  shadow  zones  behind 
ocean  obstacles,  like  seamounts,  will  be  filled  in  more  rapidly  than  by  simple  diffraction. 
Second,  this  type  of  scattering  can  lead  to  significant  transfer  of  acoustical  energy  from 
high  angles  to  low  angles.  Hence,  an  off-axis  source,  like  those  responsible  for  the  ocean 
ambient  sound  fields,  will  have  their  energy  redistributed  vertically,  thus  altering  the 
vertical  directionality  of  the  ambient  sound.  Thirdly,  scattering  along  the  front  can 
contribute  to  the  vertical  extension  of  acoustical  caustics  zones  that  abut  shadow  zones. 
Observations  reported  by  Dushaw  et  al.  (1999)  of  the  Acoustic  Thermometry  of  Ocean 
Climate  group  on  Navy  SOSUS  stations  showed  that  wavefront  data  obtained  from  these 
stations  had  experienced  dramatic  vertical  extensions  of  wavefront  caustics  several 
hundreds  of  meters  above  the  receivers.  The  resultant  scattering  caused  interference 
patterns  along  the  front  and  vertically  extended  the  front  into  shadow  zones  (Figure  1). 

Finally,  since  this  scattering  mechanism  leads  to  acoustical  energy  redistribution 
along  the  front,  the  scattering  processes  can  be  considered  to  be  coherent  in  some  regards. 
Thus,  the  scattered  energy  can  be  used  effectively  for  signal  and  array  processing  needs, 
but  the  statistics  from  this  type  of  scattering  has  not  yet  been  examined  in  detail. 

With  this  background  information  in  mind,  we  consider  a  pulse  emitted  from  a 
point  source  in  the  ocean  sound  channel.  In  this  case,  an  outward  propagating  wavefront 
is  formed.  After  traveling  approximately  1,000  km,  the  wavefront  takes  on  an  accordion 
shape  (Figure  2).  The  wavefront  is  composed  of  individual  segments  that  are  joined  by 
cusps  and  each  point  on  the  wavefront  is  identified  with  a  ray  with  a  specific  launch  angle 
of  the  source.  The  wavefront  is  thus  defined  by  the  end  points  of  all  the  rays  at  a  specific 
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time  after  emission  of  the  pulse.  After  adding  small-scale  heterogeneities  due  to  internal 
waves,  each  ray  will  be  displaced  some  distance  from  its  unperturbed  end  point.  The 
deviation  of  this  ray  end  point  can  be  either  along  the  wavefront  or  perpendicular  to  the 
wavefront.  Numerical  simulations  in  this  thesis  and  in  other  work  will  show  that  the 
deviation  is  much  larger  along  the  wavefront. 

The  primary  purpose  of  this  thesis  is  to  study  along  the  front  scattering  using  a 
hierarchy  of  models  from  a  simple  1-D,  discrete  model  to  a  fully  random  ocean  model. 
The  main  acoustical  observable  is  the  geometry  of  the  acoustical  wavefront,  and  Ray 
methods  are  used  to  quantify  the  re-distribution  of  ray  energy  along  and  across  the  front. 
Of  principal  importance  is  the  ability  to  predict  the  ocean  and  acoustic  variables  that 
could  most  affect  this  type  of  scattering.  For  example,  many  of  our  simple  reduced 
physics  models  predict  that  variance  of  along  wavefront  ray  deflections  scale  like  the 
cube  of  the  propagation  range,  while  the  variance  of  across  the  front  ray  deflections  scale 
linearly  with  range.  In  addition,  scattering  in  the  wavefront  direction  is  controlled  by  the 
vertical  gradients  of  sound  speed  while  the  scattering  in  the  ray  direction  depends  on  the 
sound  speed  perturbation  itself. 
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II.  METHODOLOGY 


Ray-based  acoustic  propagation  models  with  various  ocean  sound  speed 
fluctuation  models  were  used  to  estimate  scattering  along  and  across  the  wavefront.  The 
simplest  model  considered  was  a  discrete  1-D  model  with  small  angle  forward  scattering. 
Here  Snell’s  law  was  linearized  and  used  to  analytically  predict  the  trajectory  of  ray  paths 
through  a  series  of  random  sound  speed  layers.  The  validity  of  the  linear  approximations 
was  examined  to  determine  the  relative  contributions  of  along  vs.  across  wavefront 
scattering  as  a  function  of  the  sound  speed  variance,  the  number  of  layers,  the  width  of 
the  layers,  and  the  wavefront  angle  of  incidence. 

The  next  level  of  complexity  considered  a  continuous  2-D  random  field  of  sound 
speed  perturbations  obeying  the  (GM)  internal  wave  spectrum.  Various  analytic 
approximations  to  the  scattering  were  checked  through  the  use  of  Monte-Carlo  numerical 
ray  simulations  that  examined  the  scattering  physics.  The  simplest  case  modeled  was  ray 
propagation  through  an  internal  wave  field  that  was  anisotropic  in  depth  and  the 
horizontal  direction,  but  homogeneous  in  depth.  A  comparison  of  ray  simulation 
statistics  to  some  simple  analytic  expressions  was  easily  obtained  because  of  the  simple 
geometry. 

The  last  and  most  complicated  case  involved  a  true  ocean  environment,  with  a 
sound  channel,  and  an  anisotropic  and  inhomogeneous  internal  wave  field.  Again, 
Monte-Carlo  ray  simulations  were  performed  to  test  more  complicated  analytic 
expressions  for  scattering  along  and  across  the  front. 

A.  SOUND  SPEED  ENVIRONMENT 

In  this  numerical  experiment,  the  following  two  dimensional  sound  speed  field  is 
considered: 

c(x,z)  =  c(z)  +  Sc(x,z) ,  8c □  c(z)  (2.1) 
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where  Sc= 


f  3c  ^ 
\dzy 


C,m  and  c(z)  represent  the  sound  speed  perturbation  and  mean  sound 

f  dc^ 


speed  profile  respectively.  C,  is  the  internal  wave  displacement  field  and 


is  the 


potential  sound  speed  gradient.  The  mean  sound  speed  profile  can  then  be  modeled  in 
one  of  two  ways.  Either,  it  is  represented  by  the  constant  co=1500  ms'1  or  by  the  more 
complex  Munk  canonical  form  (Munk,  1974). 


:(z)  =  co  1+s^e 


"2(z-za  )/B 


+  2 


(z-Za)/B-lj 


(2.2) 


Specifically,  the  axial  sound  speed  (co)  is  equal  to  1500  ms'1,  the  total  water  depth 
D=5000  m,  B=1000  m,  the  sound  channel  axial  depth  (za)  was  set  to  1000  m,  and 
8=0.005515. 


First,  the  potential  sound-speed  gradient  profile  is  expressed  as 
Sc(r,t 


where  Co  is  1500  ms'1,  N  is  the  buoyancy  frequency  and  a  is  a  dimensionless  constant 
related  to  salinity  and  temperature,  approximately  equal  to  24.5. 


^ a ^ 
Sy 


N2C 


(2.3) 


From  Flatte  et  al  (1979),  a  fractional  sound  speed  variance  profile  was  imposed  in  order 
to  understand  the  internal  wave  induced  sound-speed  fluctuations: 


fi2(z)) 


=  F 


N3  (z) 


N 


ref 


(2.4) 


where  ^p‘^=6.26  x  10'8,  and  Nref  =  3  cph.  However,  near  the  surface  and  bottom,  the 
profile  of  ^p2  (z)^  will  be  modified  in  order  to  accommodate  the  internal  wave  vertical 

mode  boundary  condition  of  zero  displacement  at  those  areas.  Equation  2.4  is  important 
because  it  reveals  the  inhomogeneities  of  internal  wave  induced  sound-speed  fluctuations 
with  respect  to  depth.  As  a  result  of  a  larger  buoyancy  frequency  in  the  upper  ocean  as 
opposed  to  a  deeper  one,  these  sound-speed  fluctuations  are  larger  in  the  upper  regime. 
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The  task  now  falls  to  generating  the  random  realizations  of  sound  speed 
perturbations  that  are  induced  by  internal  waves,  where,  internal  waves  are  oscillatory 
geophysical  wave  motions  resulting  from  the  Coriolis  and  buoyancy  forces.  The  change 
in  sound-speed  due  to  internal  waves  is  on  the  order  of  1  ms'1,  with  time  scales  lying 
within  the  regime  of  the  Coriolis  and  buoyancy  periods.  The  length  scales  of  internal 
waves  are  of  order  100  m  in  the  vertical  and  10  km  in  the  horizontal.  Further,  sound- 
speed  fluctuations  caused  by  internal  waves  are  the  dominant  source  of  variability  of 
acoustic  fields  in  the  ocean. 


The  Garrett-Munk  internal  wave  spectrum  is  utilized  using  the  Colosi  and  Brown 
method  (1998).  Internal  waves  can  be  represented  as  horizontally  propagating  vertical 
modes.  Consequently,  the  internal  wave  displacement  field  C,  can  be  considered  as  the 
sum  over  all  possible  internal  waves  with  vertical  mode  number  j  and  wave  numbers  kx 
and  ky. 


Ciw  (x,y,z)  =  Re 


CO  j 

-  Jmax 

J  dkx£g(j,kx)AJ(z) 

-co  j-l 


i(kxx-CTjt) 


(2.5) 


where  Aj  is  the  jth  mode  function,  kx  is  the  horizontal  wave  number  in  the  propagation 

plane.  The  amplitudes  of  the  internal  waves  are  complex  random  variables  with  statistics 
described  by  the  GM  spectrum.  The  buoyancy  frequency  profile  is  modeled  using  a 
canonical  profile  as  shown  below. 

N(z)  =  Nse'z/B  (2.6) 


Here,  B  is  still  equal  to  1,000  m  and  the  surface  buoyancy  frequency  (Ns)  is  set  to  5  cph. 
However,  when  the  background  sound  speed  is  a  constant,  Ns=lcph  is  employed.  The 
sum  of  the  contributions  from  internal  wave  mode  numbers  up  to  jmax  =  200  was  chosen 
and  internal  waves  with  horizontal  wave  scales  from  0.5  to  409.6  km  were  utilized  in  the 
generation  of  the  random  realizations  of  internal  wave  induced  sound  speed 
perturbations. 

A  random  realization  of  the  spectrum  of  p  is  achieved  by  expressing  it  in  terms 
of  a  mode  number  and  horizontal  wave  numbers  kx  and  ky.  The  final  result  is 
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(2.7) 


G,(j\kx)  =  ^ 


N3  (z)  1  1  2 


N3  M  j2+j*  n1 


k. 


j 


-  +  - 


1 


k-+k-  2 (k;+kr) 


In 


(ki+krp+kj 

(k^-k, 


f^j 


where  k  =  — —  .  Internal  waves  are  isotropic  in  the  horizontal  direction,  but  anisotropic 

J  n0b 

in  the  vertical. 


B. 


RAY  THEORY 


The  evolution  of  the  wavefront  is  well  described  by  the  two-dimensional  ray 
equations  (Brown  et  al,  2003,  and  Beron-Vera  et  al  2003).  They  are  given  by 

(fp  1  3c 


dx  (.  2  W/2  dz 

=  c(l-c-p  )  , 


dt 


dX. 


c2p. 


dX  c  dz 


(2.8) 


where  p=sin6,/c  is  the  ray  slowness  in  the  z-direction.  The  ray  equations  were  integrated 
forward  in  time  to  Tmax=1000  km/co  =  666.667(s),  with  the  intermediate  times  being 
saved  in  increments  of  0.667  seconds.  These  equations  were  solved  numerically  using 
MATLAB. 

C.  NO  WAVEGUIDE  CASE 

For  the  case  of  constant  background  sound  speed,  simulations  were  carried  out 
with  a  source  depth  of  1,000  m  and  with  an  initial  ray  angle  of  zero  degrees.  A  total  of 
12,000  realizations  of  ray  propagation  were  computed  in  order  to  accurately  calculate  the 
necessary  statistics.  In  this  case,  a  periodic  boundary  condition  was  implemented  in 
order  to  compensate  for  rapid  scattering  to  the  ocean  boundaries.  Figure  3  shows  the  ray 
endpoints  from  all  the  realizations  at  the  following  nominal  ranges:  340,  680,  and  1,000 
km.  A  dramatic  depth  deviation  (wavefront  direction)  of  the  end  point  relative  to  the 
range  deviation  (ray  direction)  is  clearly  evident.  This  figure  dramatically  shows  the 
dominance  of  scattering  along  the  front.  The  rays  are  also  seen  to  be  scattered  into  an  arc 
instead  of  a  straight  line.  The  arc  is  fonned  because  the  rays  with  the  largest  depth 
deviation  could  not  get  down  range  as  far  as  the  rays  that  stayed  closer  to  the  original 
depth.  In  Figure  3,  the  green  curve  represents  the  parabolic  fit  to  ray  end  point 
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distributions  and  the  red  curve  shows  the  location  of  the  unperturbed  point  source  front 
(spherical  wave).  The  parabolic  fit  is  used  as  a  mean  wavefront  and  thus  this  function  is 
subtracted  off  the  ray  endpoints  to  yield  statistics  in  the  wavefront  and  ray  directions.  It 
is  clear  that  at  1000  km  the  wavefront  direction  deviations  are  on  the  order  of  tens  of 
kilometers  as  opposed  to  the  order  of  10  meters  of  the  ray  direction  deviations.  Statistical 
results  will  be  given  in  later  sections. 

D.  WAVEGUIDE  CASE 

Implementing  a  source  depth  of  1,000  m  and  an  initial  ray  launch  angle 
of6*r(0)  = +7°,  where  the  angle  Or  is  with  respect  to  the  x-axis,  the  case  of  the  Munk 

background  profile  with  internal  wave  induced  sound  speed  variations  was  simulated.  A 
total  of  12,000  realizations  of  ray  propagation  were  modeled.  Figure  4  shows  the  ray 
endpoints  from  all  the  realizations  at  nominal  ranges  of  340,  680,  and  1000  km.  The 
location  of  the  unperturbed  +7°  ray  is  denoted  by  the  green  dots.  Further,  the 
unperturbed  front  composed  of  initial  angles  between  +0°  and  + 14°  is  denoted  with  a  red 
line.  At  a  range  of  1000  km,  deviations  of  the  ray  end  points  in  the  ray  direction  are  of 
the  order  10  m.  Again,  in  the  wavefront  direction,  the  order  is  in  kilometers, 
demonstrating  the  dominance  of  along  wavefront  deviations.  Also  note  that  the 
scattering  in  the  wavefront  direction  is  seen  to  be  asymmetrical  in  the  sense  that  the 
wavefront  is  extended  more  towards  the  longer  and  not  the  shorter  ranges.  Statistics  will 
be  discussed  in  later  sections. 
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III.  1-D  DISCRETE  CASE 


The  oceanographic  wave  propagation  physics  leading  to  distinct  anisotropy  in 
wavefront  distortions  is  geometrically  complex.  To  alleviate  this  complexity,  and  better 
understand  the  precise  geometry  involved  with  this  problem,  a  plane  1-D  model  is  first 
considered.  Although,  the  model  is  not  realistic,  its  simplicity  allows  for  the  analytical 
solutions  of  the  trajectories  of  individual  rays  which  are  necessary  in  elucidating  the 
aforementioned  geometry. 

A  model  consisting  of  a  layered  random  sound  speed  medium  with  small  angle 
forward  scattering  was  considered  (see  Figure  5).  There  is  an  incident  wave  with  angle 
do  (with  respect  to  the  horizontal),  in  a  region  of  sound  speed  Co.  The  incident  wave 
refracts  via  Snell’s  Law  with  angular  deviation  567  into  region  1  of  a  slab  with  specified 
width  L,  where  Colosi  and  Flatte  (2007)  defined  the  sound  speed  asc^c^+ScpSc,  Gjj  c0. 

The  scattered  ray  then  moves  a  distance  xi  in  the  horizontal,  taking  time  fi  to  cross  the 
slab  in  this  region.  The  process  continues  in  this  manner  for  a  specified  number  of  N 
slabs. 


cos  0 

Colosi  and  Flatte  (2007)  used  Snell’s  Law, - =  constant,  to  express  a  small 


angular  deviation  when  there  is  a  small  sound  speed  perturbation  5c  ■ .  Their  results  are 
given  by 
-fij 


567  = 


tan  6,, 


(3.1) 


where  p.  =  Sc]  /  c0  and  67  =  d()  +  <567 .  After  the  scattered  ray  crosses  each  layer,  an  x- 
deflection  and  a  travel  time  result.  Their  resulting  equations  are  exact  relations. 


L 


Xj  = 


tan(6,o  +  56,j) 


tj  = 


L 


Cj  sin(6,o  +  56lj) 


(3.2) 
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After  expansion  in  small  angles,  Flatte  and  Colosi  (2007)  arrived  at  the  following 
approximation: 


Sxj  =  Xj-X0  : 


xo- 


W 


sin"  0o 


,  8tj  =  tj-to  =  tojUj 


008  26*0 
sin2  0o 


(3.3) 


Here,  the  unperturbed  deflections  of  the  ray  path  are  represented  by  the  values 
xo=L/tan6*o  andt0=L/(c0sin#0) .  According  to  the  results  (Equation  3.3),  both  the  x- 
deflection  8xj  and  the  sound  speed  perturbation  p,  will  always  have  the  same  sign.  That 
is  to  say,  should  pj  be  positive  then  80;,  must  be  negative,  which  results  in  the  ray 
having  a  larger  glancing  angle  and  thus  a  larger  Sx^.  The  opposite  also  holds  true. 

However,  the  sign  of  the  time  shift  is  dependent  upon  both  the  angle  of  incidence  and  the 
sound  speed  perturbation.  For  the  case  of  interest  0o  <  45°,  the  time  shift  Stj  always 
holds  the  same  sign  as  the  x-deflection.  With  p,  >  0,  the  path  length  is  longer  and  the 
time  shift  positive.  Thus,  the  sound  speed  increase  is  unable  to  compensate  for  the 
increase  in  path  length.  The  result  is  an  increase  in  time  for  the  ray  to  cross  the  slab.  The 
opposite  occurs  when  p  <  0,  which  means  that  the  ray  crosses  the  slab  in  less  time.  That 
is,  for  0o  <  45°,  changes  in  the  ray  path  length  dominate  the  scattering  of  the  rays. 


A  random  walk  model  was  tested.  First,  the  sound  speed  fluctuations  pn  are 
assumed  random,  within)  =  0  and(pnpm)  =  ^p2^Snm.  Next,  Colosi  and  Flatte  (2007) 
derived  the  following: 


Xn 


dx=V  8x  = 

£—t  J  n  i—t 

j=l  sin  p>  j=i 


V  cos6*0  ^ 

2 rwL2^j 

sin  0, 


(3.4) 


o  j=i 


ti  ‘  sin  6),  ~i‘ 


cos  20o 
sin3  0Q 


Em 

0  j= 


(3.5) 


where  dx  and  dt  are  the  summations  ofSxj  andStj ,  respectively.  Colosi  and  Flatte  (2007) 
derived  the  results  for  the  moments  of  dx  and  dt. 


\- 

6  8xo  6 

/- 

vsin2  0o y 

N  01- 


cos2  00 
sin6  0o 


(3.6) 
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(3.7) 


(dt2)  =  5to2 


^  cos  20o  y 
v  sin2  0o  j 


cos2  20o  L2 
sin6  0o  co2 


The  growth  of^dx2  ^  and^dt2^  is  diffusive  because  they  grow  linearly  with  the  number  of 
scattering  events  (N). 

The  approximations  for  Sx,  and  Stj  were  numerically  tested  with  random 
realizations  and  sound  speeds  p .  The  following  parameter  values  were  chosen:  sound 
speed  Co  =  1500  ms'1,  depth  of  the  layer  L  =  100,  number  of  layers  N  =  100  and 

(i 8c2^2=  0.15  ms'1.  The  variance  of  dx  and  dt  were  computed  for  an  initial  range  of 

angles  between  2  and  88  degrees.  The  approximations  proved  to  be  excellent,  as 
demonstrated  in  Figure  6. 

A.  DECOMPOSITION  ALONG  AND  ACROSS  THE  FRONT 

The  scattered  rays  need  to  be  compared  to  the  unperturbed  wavefront  by 
correcting  the  scattered  ray  to  the  same  time  as  the  unperturbed  wavefront.  With  the 
assumption  that  the  scattered  ray  propagates  at  the  angle  0o,  Colosi  and  Flatte  (2007) 
arrived  at  the  following  geometrical  analysis  (see  Figure  7): 

Aw=dx  sin0o ,  and  AR=dx  cos0o-dtco  (3.8) 

where  Aw  signifies  the  deviation  along  the  wavefront  and  Ar  is  the  deviation  in  the  ray 
direction,  which  is  normal  to  the  wavefront.  The  second  term  in  Ar  represents  the  time 
correction  back  to  the  time  of  the  unperturbed  front.  Through  Equation  3.8,  a  critical 
understanding  of  the  relative  sizes  of  the  two  wavefront  deviations  is  achieved.  As  a 
result  of  the  fact  that  dx  and  dt  have  the  same  sign  for  0O  <  45° ,  the  two  contributions  to 

Ar  nearly  cancel  and  the  rays  fall  back  along  the  unperturbed  front.  As  a  result  of  this 
near  cancellation,  Ar  is  reduced  by  a  factor  of  tan  0  relative  to  Aw .  This  is  a  significant 
result  as  tan  0  can  be  quite  small  for  launch  angles  near  zero  degrees.  Glancing 
incidences  created  a  larger  distance  and  thus  longer  time  for  crossing  the  slab, 
accentuating  the  effect. 
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The  second  moments  were  then  obtained  using  the  random  walk  model  described 


earlier. 


(a2w)  =  nl!(h!) 


COS2  60 
sin4  Oo  ’ 


(3-9) 


When  Oo  approaches  90  degrees,  or  near  nonnal  incidence,  cos2  Oo  approaches  zero, 
thus  ^ Aw2  ^  becoming  very  small;  however,  sin2  Oo  approaches  one,  thus  ^  Air  ^ 

approaches  a  constant.  As  Oo  approaches  zero  or  glancing  incidence,  ^  Aw2  \  grows  more 
rapidly  than  (Ak2  ^  because  sin4  Go  grows  more  rapidly  than  sin2  Oo .  At  <9o=45  degrees, 
the  curves  intersect  (see  Figure  8).  The  figure  also  shows  that  the  Monte-Carlo 
simulation  does  very  well  in  approximating  ^Ar  2  ^  and  (Aw2^)  and  is  very  close  to 
theoretical  estimates. 
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IV.  2-D  CONTINUUM  MODELS:  NO  WAVEGUIDE 


The  discrete  1-D  model  used  in  the  previous  chapter  did  well  to  describe  the 
geometry  of  the  problem  in  question,  but,  nonetheless,  brought  forth  serious 
disadvantages.  In  the  1-D  model,  the  x-component  of  the  ray  slowness  (cos  67c)  was  a 

conserved  quantity,  thus  it  can  strongly  affect  the  fluctuations  in  the  system.  From 
Beron-Vera  et  al.  (2003)  and  Beron-Vera  and  Brown  (2004),  in  ray  systems  in  which  the 
sound  speed  is  a  function  of  only  one  variable.  The  wave  number  component 
perpendicular  to  the  sound  speed  gradient  is  constant.  When  the  sound  speed  varies  in 
multiple  directions  the  wave  number  components  are  not  conserved  and  thus,  ray  chaos 
can  occur  (i.e.,  exponential  sensitivity  to  initial  conditions).  Further,  linearization  of  the 
equations  was  prohibited  in  the  1-D  model  at  an  incident  angle  of  zero  as  total  reflection 
issues  occurred. 

A.  DEVIATIONS  IN  THE  WAVEFRONT  DIRECTION 

In  the  case  of  no  sound  channel,  the  geometry  represents  an  unperturbed  ray  as  a 
straight  line  from  the  source.  The  direction  x  is  along  the  unperturbed  ray  starting  at  the 
point  (0,0).  If  there  are  no  fluctuations  in  sound-speed,  the  ray  arrives  at  position  z  =  0 
and  x  =  R  after  time  T,  as  opposed  to  Z  and  R+X  after  time  T  with  fluctuations  (  see 
Figure  9).  The  deviation  of  Z  and  X  is  along  the  wavefront  and  perpendicular  to  it, 
respectively.  Here,  Flatte  and  Colosi  (2007)  considered  an  unperturbed  ray  angle  at  zero, 
thus  9  will  be  used  to  symbolize  the  angle  at  position  X  caused  by  the  perturbations. 
The  perturbation  p  is  responsible  for  the  deviation  in  the  ray  angle  and  can  be  estimated 
from  the  ray  equations.  In  this  way,  Colosi  and  Flatte  (2007)  calculated  the  deviation  in 
angle  after  traveling  a  distance  dx  asd#[  -p<7x,  where  p  is  the  random  vertical 
gradient  of  fractional  sound  speed.  Their  equation  for  the  ray  angle  at  Range  R  is  then: 

^R)  =  |oR^  =  -JoRdxp(x)  (4.4) 

The  integration  path  of  Equation  4.4  is  along  the  path  of  the  unperturbed  ray, 

alongside  the  x-axis.  Thus,  after  a  range  R,  the  variance  in  the  ray  angle  is: 
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(4.5) 


(°2 ) = r  dxii0R  dx2  m »  h)  d  r dx  (^2 ) 


Lh 

i T 


/  2  \  LhR 

WlT 

Z 


In  evaluating  this  integral,  Colosi  and  Flatte  (2007)  made  the  assumption  that  the 
correlation  length  of  the  sound  speed  gradient  is  short  in  comparison  to  the  range  R  in 
equation  4.5.  The  correlation  lengths  of  the  internal  waves  along  and  traverse  to  the  ray 
are  defined  as  LH  and  Lz  respectively.  Equation  4.5  shows  that  the  rms  angle  with  respect 
to  the  horizontal  grows  as  the  square  root  of  range. 

Now  a  derivation  for  the  result  for  the  along  wavefront  scatter  is  made.  Using  the 
supposition  that  a  ray  propagates  to  some  position  xi  and  undergoes  a  small  angular 
deflection6,(xi) ,  the  scattered  ray  will  then  propagate  to  a  range  R  following  a  straight 

line  and  to  a  depth  deflection  of  dz(R,xi) .  Simple  geometry  (Figure  10)  shows  that 
,  ,  z(xi)  dz  /  \ 

tan#(xi)  □  #(xi)  =  — - - — - ,  or  — (xi)  =  R-xi  (4.6) 

V  V  ’  R-xi  ’ 


Next  Flatte  and  Colosi  (2007)  took  the  sum  of  all  deflections  so,  the  total  z-deviation  at 
range  R  is 

Z(RW„  —  (xi)rf^  =  -Jo  (R-x,)TA,.  (4.7) 


The  variance  of  Z  is  then 


i*R  pR  | 

(  Xl^ 

( 

X2  ^ 

axi  dx 2 

1- 

«_ 

o 

«— 

o 

l  rJ 

V 

Ry 

fl(xi)p  (x2)\. 


(4.8) 

The  two  integration  variables  are  within  Lh  of  each  other  due  to  a  restriction  imposed  by 
the  correlation  function  ofp  .  Thus,  Colosi  and  Flatte  (2007)  reduced  the  integral  to  the 


following  result: 

Z2(R))DR2J> 


Xl^2 

'r  j 


(Rf) 


Lh  □ 


R3(p2)Lh 

3U 


(4.9) 


As  is  evident,  nns  deviation  along  the  wavefront  grows  as  the  3/2  power  of  the  range. 
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Numerical  evaluations  can  now  be  perfonned  for  Qnm  andZrms  through  the 
estimation  of  the  internal  wave  factors  of  Equations  4.5  and  4.9.  The  following  rms 
values,  6 rms  3, 6  and  Zi®  □  60  km  are  a  result  of  propagation  to  a  range  of  1000  km  with 
sound  speed  fluctuation  values  characteristic  of  internal  waves:  Lh  =  10  km,  Lz  =  0. 1  km, 
and  (d2)  =  10~8 . 


B.  DEVIATIONS  IN  THE  RAY  DIRECTION 

The  perturbed  time  5T  can  be  calculated  along  the  ray.  X  is  given  bycoST, 
where  8T  is  the  travel-time  deviation  for  arrival  at  R  calculated  along  the  unperturbed 
ray.  Colosi  and  Flatte  (2007)  arrived  at  the  following: 

z  \  pR  co  r R  /  \ 

X=coT(R)-R=jo  — - -dx-RD  jo  p(x,z)dx  (4.10) 

c  ( X,z ) 

The  variance  of  X  is  then: 

(x2)  =  (p2)RLH  (4.11) 


Comparing  rms  values,  deviation  along  the  wavefront  (Z  y2  has  an  rms  value  of 

i 

60  km  in  contrast  to  an  rms  value  of  only  (x2^  =  10  m  of  the  deviation  along  the  ray. 


This  is  a  surprising  result.  It  says  that 


1/2 


X" 


1/2 


□  1000,  stating  that  deviation  is  one- 


thousand  times  greater  along  the  front  than  across  it. 


C.  EXPRESSIONS  FOR  SOUND  PROPAGATION  THROUGH  INTERNAL 
WAVES 

The  wavefront  wanderr2  (Flatte  et  ah,  1979),  is  directly  related  to  the  deviation  of 
the  wavefront  along  the  ray  (x2  'j .  As  such,  Colosi  and  Flatte  (2007)  derived  the 
following: 
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(4.12) 


(x2)  =  c2r2  =  JoR(h2(z))lp(^i=0,z)Jx 

where  Lp  is  the  correlation  length  of  the  internal  waves  in  the  x-direction  (along  the  ray) 
and  6*.  is  the  ray  angle  with  respect  to  the  horizontal.  Also,  Lp  is  dependent  upon  6r 
because  internal  waves  are  anisotropic.  The  deviation  in  ray  angle#  is  different  in  that  it 
is  not  dependent  upon  p ,  but  instead  upon  the  correlation  of  the  vertical  gradient  of  p . 
Therefore, 

(s!WoV(z))L'p(ft=o'z)‘fc  (4-B) 

where  Lp  represents  the  correlation  length  of  the  vertical  gradient  of  the  internal  waves  in 

the  x-direction.  Appendix  A  gives  expressions  for  these  correlation  lengths.  Colosi  and 
Flatte  (2007)  also  show 

(z2)=R3f  (‘-fj  (^(z)K(^°-z)A  <4-14) 

D.  NUMERICAL  RESULTS 

Figure  1 1  shows  the  statistics  of  the  nns  ray  angle  and  of  the  rms  along  the 

\_  _1_ 

wavefront  (z 2^2  and  across  wavefront  (x2^2  deviations.  With  regards  to  the  range 
scaling,  the  rms  along  the  front,  plotted  in  the  upper  portion  of  Figure  10,  behaves  as 

i 

expected,  rapidly  growing  as  R3/2.  Similarly,  the  rms  ray  angle (#2)2  displayed  in  the 
lower  portion  of  the  plot  grows  as  expected  by  R  .  However,  the  across  the  wavefront 
statistics  ^X2^2  show  some  deviation  from  the  prediction  of  Equation  3.8.  Despite  the 

subtraction  of  the  parabolic  curve  from  the  Monte-Carlo  simulation,  at  ranges  beyond 

1 

/  2  \ —  •  1  /2  •  3  /  2 

300  km,  the  growth  rate  of  (X  V  ceases  to  increase  by  R  and  jumps  to  a  rate  of  R  for 

\_ 

reasons  not  known  (future  research).  Consequently,  a  new  approach  for  (x2^2  is 
required  for  R>300  km. 
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E. 


OCEAN  APPLICATIONS:  HORIZONTAL  PLANE  SCATTERING 


An  important  result  is  obtained  when  an  examination  of  the  horizontal  spreading 
of  acoustic  rays  in  the  horizontal  (lat,  long)  plane  of  the  previous  section  is  performed. 
The  effect  is  that  this  out  of  plane  scattering  may  now  play  an  important  role  should  the 
wavefront  become  obstructed  by  an  island  or  seamount,  which  would  lead  to  the  healing 
of  the  wavefront. 


Colosi  and  Flatte  (2007)  modified  Equation  4.9  to  yield  the  horizontal  deviation 
Y2 )  for  an  isotropic  sound  speed  fluctuation  field: 


(4.15) 


An  order  of  magnitude  can  now  be  obtained  using  parameters  (same  as  the 
previous  section)  typical  of  ocean  internal  waves  and  a  propagation  range  of  1000  km. 

The  rms  value  of  the  out  of  plane  deviation  (y2^2  □  600  m,  which  shows  the  effects  of 
internal  waves  to  appear  quite  small.  However,  for  a  random  field  of  mesoscale  eddies 
with  10  5  and  Lh  □  50  km,  a  much  larger  nns  deviation  of  10  km  is  obtained. 

Diffraction  grows  as  R  as  opposed  to  scattering  healing  which  grows  as  R  .  That  is, 
diffraction  dominates  in  shorter  ranges  and  scattering  healing  dominates  in  larger  ranges. 


F.  WAVEFRONT  HEALING  IN  THE  VERTICAL 

The  following  images  show  the  progression  of  the  healing  of  a  wavefront  in  the 

vertical  past  an  underwater  obstruction  with  respect  to^Z2^2 .  A  sound  source  emitting 

rays  is  considered.  Some  of  the  rays  will  hit  the  obstruction  and  others  will  go  over  it. 
As  a  result,  a  void  or  shadow  zone  will  occur  with  respect  to  the  vertical.  As  such,  the 
along  front  scattering  in  question  is  demonstrated  to  reveal  faster  infilling  of  the  shadow 
zones  than  predicted  by  regular  diffraction  (not  shown).  At  only  100  km  past  the 
obstruction,  the  wavefront  completely  fills  in  the  shadow  zone. 
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V.  INCORPORATION  OF  THE  WAVEGUIDE 


The  results  of  the  no  waveguide  model  were  very  encouraging,  but  now  we  want 
to  treat  the  important  effects  of  the  waveguide.  Therefore,  it  was  necessary  to  modify  the 
theory  to  incorporate  the  ocean  waveguide  and  depth  dependent  ocean  internal  waves. 
To  better  understand  this  new  case  a  geometrical  approach  is  taken.  In  the  no- waveguide 
case,  the  ray  propagates  to  some  point  x’,  scatters  and  then  propagates  to  the  final  range 
R.  From  Snell’s  law,  that  ray  will  shift  some  vertical  differential  distance  d#  after 

dz 

traveling  some  horizontal  distance  R-x’.  As  a  result,  —  equals  R-x’  (see  Figure  10). 

d  6 

The  geometry  in  the  waveguide  case  is  similar  to  the  no-waveguide  case  except  that  the 
ray  is  now  following  a  curved  ray  path  and  not  a  straight  line  (Figure  13).  Thus,  the 
derivations  for  the  variance  in  X,  Z,  and  9  remain  relatively  the  same  with  the  exception 
dz 

that —  is  now  dependent  upon  the  ray  tube  function^  .  Colosi  and  Flatte  (2007)  show 
d  6 

that 


dz  ^2(x’) 

d9y  ’  q2'(R) 


(5.1) 


The  following  statistical  equations  from  Colosi  and  Flatte  (2007)  are  for  the  waveguide 
case: 


f ds  lp  )’(»2)= JT ds  )  F  •  A  )■ =  f ds  <^2 )  k 


02(R)t 


(5.2) 


A.  NUMERICAL  RESULTS 

In  comparison  to  the  results  of  the  previous  chapter,  the  nns  ray  angle 

statistics'#2^2  continues  to  show  the  expected  growth  of  R1/2  with  some  deviation  at 
larger  ranges  followed  by  a  tapering  off  at  a  range  of  about  1000  km.  The  across  the 

wavefront  statistics  (x2^ 2  now  grow  as  R1/2  throughout  the  entire  range  interval  as 
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'I/') 

opposed  to  the  R  ~  seen  in  the  no-waveguide  case  after  300  km  range.  Theory  fits  well 

\_ 

for^Z2V  up  to  600  km,  but  shows  deviations  at  greater  range.  For  large  ranges,  there 
maybe  a  problem  with  the  calculation  of  the  ray  tube  function  C,  past  ranges  of  600  km. 

3  /2 

Whereas,  the  Monte-Carlo  simulation  tapers  off  at  larger  ranges  and  grows  nicely  at  R  , 
the  theory  shows  a  dramatic  deviation  from  the  model  statistics  (Figure  14). 
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VI.  CONCLUSIONS 


The  ray-based  physical  models  did  very  well  in  elucidating  the  physical 
mechanism  for  the  ocean  scattering  processes  that  tend  to  redistribute  acoustic  energy 
dominantly  along  the  wavefront.  The  models  in  question  do  a  relatively  good  job 
predicting  the  wavefront  distortions  caused  by  the  ocean  internal  waves  with  respect  to 
Monte-Carlo  simulations.  For  the  simple  1-D  model,  dx  and  dt  have  the  same  sign,  thus 
ray  direction  variance  is  diminished.  Further,  wavefront  scattering  is  dominated  by  path 
length  changes  of  the  ray.  In  the  continuum  models,  ray  direction  deviations  depend  on 
jll  and  grow  diffusively  as  R  .  Also,  wavefront  direction  deviations  depend  on  p  and 
there  is  a  moment  arm  effect  that  makes  (z1^  grow  like  R32.  With  regards  to  the 

waveguide  case  model,  the  problem  with  the  theoretical  prediction  of  the  Z  statistics 
showing  the  dramatic  deviation  with  the  Monte-Carlo  simulation  must  still  be  worked 
out.  Nonetheless,  the  statistics  provided  for  seven  degree  angle  rays  propagating  through 
Garrett-Munk  internal  waves  in  a  Munk  sound  speed  profile  remain  very  compelling,  and 
that  the  distribution  of  acoustic  energy  along  the  wavefront  is  dominant. 
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VII.  FUTURE  APPLICATIONS 


This  paper  looked  only  at  one  ray  with  a  launch  angle  of  seven  degrees  using  a 
ray  theory  approach.  This  method  has  some  obvious  inherent  problems.  First,  a  source 
will  emit  more  than  one  ray;  and,  the  present  ray  theory  does  not  give  any  information 
concerning  the  amplitude  of  the  acoustical  field.  To  apply  the  results  presented  herein  to 
an  operational  environment,  the  energy  distribution  of  all  the  rays  must  be  inferred  by 
observing  where  the  ray  endpoints  end  up  at  along  the  wavefront.  Further,  the 
calculations  should  also  be  compared  with  those  implementing  parabolic  equations  and 
nonnal  mode  theory  to  verify  results  and  the  effects  of  diffraction.  This  new  approach 
will  help  to  understand  the  practical  implications  of  along  the  wavefront  scattering.  One 
application  is  the  rapid  infilling  of  shadow  zones  behind  underwater  obstructions,  the 
significant  transfer  of  acoustical  energy  from  high  to  low  angles,  and  the  dramatic 
vertical  extension  of  acoustical  caustic  zones  that  abut  shadow  zones  that  were  reported 
by  the  Acoustic  Thermometry  of  Ocean  Climate  group  at  Navy  SOSUS  stations  in  1999 
will  be  more  completely  elucidated. 

For  wavefront  healing  the  important  result  is  the  comparison  of  diffractive  healing 
which  scales  as  R  and  scattering  which  scales  as  R  “.  At  short  ranges,  diffraction 
dominates,  but  at  longer  ranges  scattering  will  dominate.  Where  this  transition  occurs  is 
unknown. 

Further,  the  theory  needs  to  be  compared  to  experimental  data.  Because  a  beam 
ensonifies  only  a  small  group  of  launch  angles,  an  experiment  involving  the  transmission 
of  a  narrow  angle  beam  using  an  array  of  sources  can  meet  this  need.  Littoral  Warfare 
Advanced  Development  (LWAD)  would  be  an  appropriate  platfonn  to  conduct  this 
experiment. 
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APPENDIX  A.  CORRELATION  LENGTHS:  THE  GARRETT- 

MUNK  SPECTRUM 


The  GM  spectrum  is  defined  by  two  key  variables:  vertical  mode  number  j  and 
frequency  a.  For  every  vertical  mode,  there  is  some  depth  they  are  dependent  upon. 
These  depths  are  in  turn  locally  represented  by  a  hannonic  function  with  a  wavenumber 
given  by: 

kz=ajN(z)  (A.l) 


a  =  ij>N(z) 


(A.2) 


where  D  is  the  ocean  depth  and  N(z)  is  the  buoyancy  frequency  profile.  The  following 
equations  are  expressions  for  the  correlation  length  of  p  in  the  direction#,  for  the  GM 
model: 


LP(#,z)  =  (  j"1  ^  — -  — - —  F  i 


n  acoi 


.  \  Jmax  l  1 


N(z)tan# 


J=1  J  J  +  J* 


Mj 


COi 


Jmax  1 

I 


x 

J 

V1 


2  ,  ;2 


j=1  J  +  J 


(A.3) 


(A.4) 


J 


/  ?  \l/2 

/  \  (2  X2  /  ?  \ -3/2  (X  +1)  +1 

,(x)  =  (x  +1)  +y(x-+l)  In  , 

7  (X"+1) 


(A.5) 


where  jmax  is  a  cutoff  mode  number  and  or  is  the  inertial  frequency,  ^li2  'j  and  N,  the 
buoyancy  frequency,  are  both  functions  of  depth  and  j*  is  usually  equal  to  three. 

The  correlation  functions  for  the  vertical  derivative  of  p  in  the 
direction# are  obtained  by  giving 

f  N(z)#a 


Up(0,z)=(aN)I(j>l— F, 


n  acoi 


v  COi  y 


(A.6) 
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1 


(A. 7) 


The  cutoff  mode  number  jmax  is  of  importance  forLp  given  the  divergence  of  (j) 
for  jmax  — »  oo  .  Through  observation,  jmax  is  established  when  the  vertical  wavenumber  of 
the  internal  waves  equals  0.1  cpm  (Duda  et  al  1992). 
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APPENDIX  B  FIGURES 
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Figure  1 .  Wavefront  data  collected  from  Navy  SOSUS  receivers. 
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Range  (km) 

Figure  2.  A  ray  computation  showing  an  acoustic  wavefront  at  roughly  1,000-km  range, 
resulting  from  a  point  source  on  the  sound  channel  axis.  The  calculation  was 
carried  out  using  a  Munk  canonical  sound-speed  profile  without  internal- wave 

sound-speed  perturbations. 


28 


Figure  3.  Range  evolution  of  the  ray  endpoints  for  12,000  realizations  of  ray 

propagation  through  internal  wave  induced  sound  speed  perturbations.  For  all 
realizations  the  initial  angle  is  zero  degrees  and  the  initial  depth  is  1 ,000  m.  The 
solid  vertical  lines  show  the  location  of  the  unperturbed  wavefront,  assumed  a 
plane  wave.  The  dash  red  curve  shows  a  parabolic  fit  to  the  ray  endpoint 
distributions,  and  the  green  curve  shows  the  location  of  the  unperturbed  point 

source  front. 
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Figure  4.  Range  evolution  of  the  ray  endpoints  for  12,000  realizations  of  ray 

propagation  through  internal  wave  induced  sound  speed  perturbations,  with  the 
Munk  canonical  profile.  For  all  realizations  the  initial  angle  is  7  degrees  and  the 
initial  depth  is  1000-m.  The  red  curve  shows  the  unperturbed  wavefront  from  a 
point  source,  and  the  green  dot  shows  the  location  of  the  unperturbed  7  degree 

ray. 
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Figure  6.  Comparison  of  statistics  for  dx  and  dt  as  a  function  of  the  initial  angle  0o , 
using  the  random  walk  model.  The  blue  curves  represent  the  Monte  Carlo 
simulations  and  the  red  curves  are  the  theoretical  estimations. 
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Figure  7.  Schematic  of  the  geometry  of  the  projection  of  the  scattered  ray  onto  the 

unperturbed  wavefront. 
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Figure  10. 


z 


Geometrical  description  of  the  application  of  Snell’s  Law  to  the  no-waveguide 

case. 


36 


<_x>  (69P) 


Figure  1 1 .  Computations  of  rms  ray  deflections  for  zero  angle  rays  propagating  through 
GM  internal  waves  (solid)  and  theoretical  expectations  (dash).  The  rms  angle  is 
diplayed  in  the  lower  portion  of  the  plot,  showing  the  expected  growth  by  R  . 
The  rms  along  the  wavefront  deflection  is  diplayed  in  the  upper  portion  of  the 
plot,  also  showing  the  expected  R  “  growth.  The  ray  direction  deflection  is 
displayed  in  the  middle  section  of  the  plot  where  the  numerical  result  has  been 

i  n 

corrected  for  the  parabolic  spreading  of  the  ray  endpoints.  The  R  scaling  breaks 

down  after  300-km  in  range. 
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Figure  12. 
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Progression  of  along  the  wavefront  scattering  past  an  underwater  obstruction 
demonstrating  a  complete  infilling  of  the  shadow  zone  at  a  range  of  100  km. 
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Figure  13.  Geometrical  description  of  the  application  of  Snell’s  Law  to  the  waveguide 

case. 
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Figure  14.  Computations  of  rms  ray  deflections  for  7  degree  angle  rays 

propagatingthrough  GM  internal  waves  in  a  Munk  sound  speed  profile.  Monte- 
Carlo  simulation  results  are  shown  with  solid  lines  while  theoretical  results  are 
dashed.  The  rms  ray  angle  is  in  the  lower  portion  of  the  plot  and  the  expected 
display  of  R1/2  is  evident.  The  rms  along  wavefront  deflection  is  displayed  in  the 
upper  portion  and  the  expected  growth  of  R3/2  is  evident.  The  rms  across  the 
wavefront  is  in  the  middle  of  the  plot  and  shows  the  expected  growth  of  R  . 
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